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Abstract. The propagation of shear Alfven wave packets in inhomogeneous magnetic fields, at the origin of their 
distortion, regardless of the occurrence of non-linear coupling, is considered. It is shown that the distortion 
mechanism can be regarded as mixing process and hence, standard "phase mixing" corresponds to the effect of an 
"Alfvenic" shear flow while enhanced dissipation at a magnetic X-point corresponds to mixing by an "Alfvenic" 
strain flow. The evolution of the global wave field is supposed to result from the dynamics of a superposition of 
wave packets and a kinetic equation for the wave energy is obtained following this eikonal (WKB) description. 
Since shear Alfven wave packets experience continuous shearing/straining while transported by an inhomogeneous 
Alfvenic flow Va, their mixing process, in physical space, is also a cascade of wave energy in k-space. The wave 
energy spectrum resulting from this linear mechanism of energy transfer is determined for the special case of 
waves propagating along chaotic magnetic field lines, the analog of a chaotic mixing process. The latter follows a 
fc _1 power-law, in the energy conserving range in k space. 
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1. Introduction 

Phase mixing, the propagation and enhanced dissipation 
of shear Alfven waves, remain active areas of research in 
studies of the solar corona, both in relation with its heat- 
ing mechanism and the acceleration of the fast solar wind; 
see e.g. reviews on the subjects (Browning 1989; Walsh 
and Ireland 2003; Ofman 2005) and references within. The 
reason is that if Alfven waves are excited and dissipate 
some of their energy in the solar corona, this contributes 
to the overall coronal energy budget. 

Due to the very low collisionality of the medium, a 
main difficulty of waves dissipation theories is, however, 
the justification of how these fluctuations can thermalizc, 
i.e. reach the very small dissipative scale, before they leave 
the corona. But since the coronal plasma is highly inho- 
mogeneous, it is quite natural to invoke waves interaction 
with the ambient background as a reason for the creation 
of fine enough scales in the wave field. It is in this con- 
text that, phase mixing of Alfven waves was proposed, by 
Heyvaerts and Priest (1983), as an enhanced dissipation 
mechanism. They considered simple inhomogeneity of the 
medium, in one direction perpendicular to the unidirec- 
tional background magnetic field. The existence of a shear 
in the Alfven velocity results in the creation of pro- 
gressively smaller transverse scales in the wave field while 
it propagates and is transported along the unperturbed 



magnetic field. In this case, the build-up of a diffusive 
scale is linear in time [dk/dt oc est), resulting in a dissipa- 
tion time proportional to S 1 ^ 3 « S. Here, S is the large 
Lundquist number S = t^/ta ^ 1 of the plasma. The 
corresponding diffusive scale is S" 1 ^ 3 <C 1. Therefore, the 
dissipation time is found to me much shorter than the 
very large diffusion time t v , proportional to S 3> 1, for 
the wave to dissipate in the absence of inhomogeneity. 

More recent studies have shown that the creation 
of a dissipative scale can be made exponentially fast 
(dk/dt oc cst.k), with associated time of the order of 
In 5 << S 1 / 3 « S (Similon and Sudan 1989; Petkaki 
et al. 1998; Malara et. al 2007). The latter improved effi- 
ciency of the dissipation results essentially from the con- 
sideration of a more complex magnetic geometry, with the 
property of local exponential separation or divergence of 
neighboring magnetic field lines. In this case, the diffusive 
length scale is proportional to S^ 1 ^ 2 . Therefore, fast dis- 
sipation is possible in three-dimensional chaotic magnetic 
fields (Similon and Sudan 1989, Petkaki et al 1998, Malara 
et. al 2000, Malara et. al 2007) and in two-dimensional reg- 
ular ones (Smith et al. 2007, Malara et al. 2003, Ruderman 
et. al 1998), provided they possess a certain degree of com- 
plexity. 

One main goal of this work is to show that the propa- 
gation and distortion of shear Alfven wave packets, in in- 
homogeneous magnetic fields, can be regarded as a "mix- 
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ing" effect. As is well known, the kinematics of mixing are 
governed by the topological properties of the flow. It nat- 
urally follows that the distortion of Alfven wave packets, 
propagating in complex magnetic fields, is the counterpart 
of a chaotic mixing process(Ottino 1988). 

The interaction among shear Alfven waves is also a 
well-studied mechanism of energy transfer to small scales, 
the non-linear distortion resulting from elastic collisions 
between counter-propagating wave packets producing a 
flux of their energy to small dissipative scales. Yet, lin- 
ear distortion through interaction with the background 
produces a similar flux. In the following, we neglect non- 
linearities in the magnetohydrodynamics equations. The 
reason is to focus primarily on the spectral properties 
of this energy transfer for waves propagating in complex 
magnetic fields. 

We obtain a kinetic equation for the wave energy from 
a representation of the global wave perturbation as an en- 
semble wave packets (Similon and Sudan 1989). Hence, 
the wave energy spectrum F(k,t) is the sum of the en- 
ergy of each individual wave packet, with associated wave 
vector and position, whose time dependence is given by 
Hamiltonian ray equations. This WKB method has al- 
ready been employed by several authors (Petkaki et al 
1998, Malara et. al 2000, Malara et. al 2007) to study 
shear Alfven waves propagation and dissipation in differ- 
ent types of coronal magnetic structures. Our study is a 
continuation of these works. In section 4, we determine the 
transient and stationary energy spectra, F(k,t) and F(k), 
for the special case of non-interacting Alfven waves propa- 
gating in a synthetic chaotic magnetic configuration with 
quasi-uniform rate of field lines exponentiation (Similon 
and Sudan 1989). In section 3, we show the correspondence 
which exists between the problem of the enhanced dis- 
sipation of Alfven waves, transported by inhomogeneous 
Alfvenic flows, and the problem of mixing (Ottino 1988), 
reviewed in section 2. Conclusions are given in the last 
section of this article. 

2. Mixing of a passive tracer 

The mixing of a passive field ip advected by a flow U is 
described by the linear advection-diffusion equation 



dip 
~dt 



U.VV> = r)V 2 4> 



(1) 



In absence of flow, U = 0, the field ip only diffuses on a 
large time scale, for small diffusion coefficient 77, because 
this time is proportional to r) . If the flow has some spa- 
tial variation, i.e. it is inhomogeneous, the advected field 
is stirred. It then develops increasing gradients, until dif- 
fusion takes over and ultimately mixes the field. Because 
diffusion alone is slow to mix, a stirring flow is required to 
accelerate the dissipative process. An important issue is 
therefore to find flows U that have the property of being 
good mixers. For the enhanced dissipation of shear- Alfven 
waves, this property concerns the Alfvenic flow V^, which 
by definition, is responsible for the advection of the wave 



field. To show this, we first concentrate on two standard 
and simple flow configurations having their counterpart in 
the problem of dissipation of Alfven waves, and compare 
their mixing efficiency. 

Consider the effect of shear flow U = (UqU, 0, 0), 



dip 1 dip d 2 ip d 2, tjj d 2 ip 

~at + ° y dx~ = v{ d^ + d^ + d^> 



(2) 



with initial condition ip(x,y,t = 0) = tpo sin kox, which 
has no dependence in y. When 77 = 0, the field is just 
advected, resulting in ip(x,y,t) — sinko(x — U^yt), i.e an 
unbounded creation of smaller scales along y, the shear- 
wise direction. When 77 7^ 0, a dissipative cut-off is intro- 
duced. In this case, we can look for a solution of the form 
ip = ipof(t) sin ko(x — U' Q yt), which gives 



/(t) = exp(-»7*g(t • 



U' 



.3 



(3) 



hence, the initial field evolves according to 

U' 2 

ip{x,y,t) =iP smk {x-U^t)eM-vk 2 o{t+-§-t 3 )) (4) 

From the above solution, we obtain that, for 77 <§; 1, the 
mixing time scale associated with the shear flow is 



7 "shear ^ {V^Q kg ) ^ 



(5) 



The associated dissipative length Id is found from Eq. (J2J) , 
by balancing the advection term, the second on its left, 
with the diffusion (along y) term, the second on its right, 
Vfaty/L) ~ V (iP/Q): 



h ~ (777) 



U' 



(6) 



By performing a formal spectral transformation of Eq. 
it is also found that individual modes evolve according to 



dipk 
dt 



U n k>, 



dipk 



-V(kl + k 2 v )ip k , 



(7) 



y dk x 

with the right-hand side of the above expression represent- 
ing advection in fc-space. This also indicates that a shear 
flow results in a spectral expulsion towards large absolute 
value shear-wise wave numbers, corresponding to an in- 
crease tilting of the phase front in the flow-wise direction. 
Hence, the physics behind the accelerated diffusion of the 
tracer ip in a shear flow is intuitively analogous to the 
one of phase mixing, which occurs for = (V^j/,0, 0). 
To gain a further insight into the problem, let us trans- 
form the advection-diffusion equation, by introducing the 
"shearing coordinates" that follows the streamlines : 

[x 1 = x-U' yt : y' = y,t' = i\ (8) 

From the chain rule for derivatives : 

d d d . d d d , d d 

[ dx~ = dx^ , dy~ = ~ Uot dx~ + W dt = ~ UoV dx~ + dF 1 ' (9) 



the advection-diffusion is now described by 



(10) 
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Keeping only the dominant terms in this equation, leads, 
after introduction of the new variable t = (f/ 2 i 3 )/3, to a 
diffusion equation : 



d 2 x 



(11) 



As is well known, there are two self-preserving solutions of 
a diffusion equation: harmonic and gaussian. Their tem- 
poral evolution is different, however, their damping time is 
identical. Indeed, contrary to the harmonic solution stud- 
ied above, an initial gaussian packet has its amplitude that 
decays as l/(^r) 1 / 2 , while it spreads. For the phase mix- 
ing problem, such power law behavior of localized initial 
wave-packets, i.e its diffusion along field lines, was noticed 
and studied by Hood et al. (2002); see also Tsiklauri et al. 
(2003). 

Now let us consider passive advection in a simple strain 
U = (—ax, ay, 0) with harmonic initial condition indepen- 
dent of z: 



dt 



(hp 
dx 



dip .d 2 ip 
aV dy- = ^fe 



d 2 j> 
d 2 y 



d 2 z 



A solution of the form 

= ipof(t) sin(k x (t)x + k y (t)y) 

is sought for, and hence, 



k x a; 



-k y a, 



k x (t) = ko expat; k y (t) = kg exp —at, 



f = -r,{kl + k 2 y )f 

/(t) = cxp[-^(exp 2 « t 
la 



-1))] 



(12) 

(13) 

(14) 
(15) 

(16) 
(17) 



The field amplitude has super-exponential decay in time. 
It follows that for a > 0, the mixing time due to a strain 
Tstram is reduced compared to T shear , i.e. 



1~str 



oc (1/a) ln(a/77fc 2 ), 



the reason is, that for a strain, the wave vector increases 
exponentially with time and not linearly. And we can an- 
ticipate that this situation corresponds to enhanced dis- 
sipation of Alfven waves in the presence of a magnetic 
X-point configuration with \a = (ax, ay,0). This result 
remains unchanged for the three dimensional stagnation 
flow U = (ax, ay, —2az). The dissipative scale length for 
advection and diffusion to balance each other in Eq. (fT2"l) . 
is given by 



Id 



(-) 1/2 
a 



(19) 



The general problem of the mixing of a passive scalar, 
in arbitrary flows, can be approached from a WKB de- 
scription (Antonsen 1996). This is done, for instance, by 
introducing a filter, ip(k, x, t) — J d 3 x exp[— x 2 /A 2 +ik(x— 
x')]?/>(x', t), with support A <C L, L being the length scale 
of variation of U. It is a local Fourier transform in x. 



Applying this transformation to the advection-diffusion 
equation fl}, we obtain 



D t ^p = —rjk 2 tp 



(20) 



with total derivative D t = <9 t + xV + kVfc and k 2 — 
k 2 + k 2 + k 2 . Equation (f20|) . which describes transport 
in (x, k) space is also known as the kinetic wave equa- 
tion in plasma physics (Sagdeev and Galeev 1969) and 
includes non-linear coupling terms on its right-hand side. 
It is equivalently obtained by use of the Wigner distri- 
bution (Bastiaans 1979, Antonsen 1996). The respective 
characteristics are : 



x = U = VfciJ 



k = -V(k.U) = -VH 



(21) 
(22) 



with H = k.U. The general solution of the kinetic wave 
equation (|20[) can be expressed in term of phase space 
trajectories £(x',i) and K(x',k',i) where d^/dt = v(£,i) 
and dn/dt = -Vv(£,t).K, and £(x',0) = x', k(x', k', 0) = 
k'. The solution is ^(x,k,i) = / d 3 x'd 3 k'$(x' , k', 0)<5(x - 
£(x', t))S(k - k(x', k', t)) exp(-r/ / * K 2 dt'). 

The physical interpretation of the WKB approach to 
mixing is clear. It means that the advected field is decom- 
posed into an ensemble of "packets" with length scale, 
A, much smaller than the one, L, characterizing the spa- 
tial variation of the flow U. Therefore, each packet ex- 
periences, along its trajectory, a velocity field U which 
is a linear function of the space coordinates :U(xp(t)) + 
dU/dx |xp(t) - x 7 x pW being the position of the packet 
at time t. An important quantity, therefore, is the local 
rate of strain tensor, <9U/<9x, of the flow. To sum-up, the 
packets are deformed only by the shear and the strain con- 
sidered below. This brings us to the analogous effect for 
Alfven wave packets. 



(18) 3. Mixing of shear Alfven waves 



The rays equations of optics are x = V^w and k = — Vw, 
with w(k) the wave dispersion relation. The eikonal de- 
scription of MHD waves was introduced by Weinberg 
(1962); see also Similon and Sudan (1989), Petkaki et al. 
(1998), who used it in the context of wave dissipation in 
the corona. For the non-dispersive shear Alfven waves, 
with uj = ±k.VA, the rays equations read 



x = ±V A , 

k = T V(k.V A ). 

with the energy equation being 

e± = -r]k 2 e±, 



(23) 
(24) 



(25) 



The latter is required to study the energetics of the mix- 
ing process. The equation of motion Eq. ([23|) states that 
shear- Alfven wave packets propagate along magnetic field 
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lines so they disperse as the field lines do. It provides in- 
formation on the precise path along which energy is de- 
posited. The enhanced dissipation time follows from equa- 
tion (|2~4"|) , giving the rate of stretching of wave packets, 
i.e. the variation of their wave-vector, when combined 
with the energy equation. Dispersion and stretching are 
however related, due to the following consideration of 
two infinitesimally separated trajectories and the phase 
gradient, Wtp = k, they generate in an inhomogeneous 
Alfvenic flow V^. Indeed, constancy of the phase along a 
trajectory implies that ip{£,{x + r,t),t) — -0(£(x, t),i) = 
5£(x,i).V^(£(x,t),t) = const with <5£(x,t) = ££(x + 
r, t) — <5£ (x, t) being the differential separation vector. 
Thus, d(5£.k)/dt = 0. According to Eq.([23]), the trajec- 
tory of a wave packet is £ = ±V A (^,t) and therefore 
d5£(xL,t)/dt = ±^.VVi(^(x,t),t). Hence, the stretching 
properties of the Alfvenic flow, represented by Eq. lpM)) for 
the wave vector are related to the dispersion of magnetic 
field lines. Considering the phase mixing situation, with 
an Alfvenic shear flow Va = V A yx, dk y /dt = k x V' A with 
solution k y {t) = k x V A t and integrating the energy equa- 
tion, the typical time scale for the phase mixing is found 
to correspond indeed to r s /j ear , i.e 

Tmix oc S 1/3 <C S (26) 

with S the Lundquist number S 3> 1. In the same way, for 
the enhanced dissipation of shear Alfven waves at a mag- 
netic X-point, which corresponds to a an Alfvenic strain 
flow Va — ctxn.— ctyy, a = V A , the component of the wave 
vector in the direction of contraction x, increases now ex- 
ponentially with time, k x (t) = fco exp(ai). Integrating the 
energy equation results in the above mentioned superex- 
ponential decay and, therefore, the enhanced dissipation 
time corresponds to T stra in, 

W ocln^S 1 / 3 (27) 

This dependance of the dissipation time with Lunquist 
number, for Alfvenic strain, is also typical of wave damp- 
ing in complicated magnetic geometries, either if the mag- 
netic field is stochastic (Similon and Sudan 1989) or reg- 
ular, but possesses separators (Malara et al. 2003; Malara 
et al. 2007). The physics is however the same, and is re- 
lated to the divergence of nearby magnetic field lines. In all 
these cases, the dissipation time of a wave packet weakly 
depends on the Lunquist number, only as a logarithm, 
Tmix oc In S, and not as a power law, T m i x cx S* 1 / 3 obtained 
for unidirectional magnetic field, i.e. standard phase mix- 
ing. It is now clear that the propagation of shear Alfven 
waves in an inhomogeneous magnetic field is equivalent to 
a mixing process. Since the kinematics of mixing are gov- 
erned by the topological properties of the flow, it results 
that the enhanced dissipation of Alfven waves propagating 
in a chaotic magnetic field is equivalent, from the energet- 
ics viewpoint, to a chaotic mixing process. The spectral 
property of the wave energy flux, in this linear process, is 
the subject of the following section. 



4. wave energy spectrum in chaotic magnetic 
fields 

Since within the eikonal formulation the global wave per- 
turbation can be treated as a superposition of wave pack- 
ets with wave- vector kj , each carrying an energy ej (t) , we 
can also write the wave energy spectrum in the following 
form: 

F(k,t) = 5>,(M) = 5>i(t)*(A;- |k,-(i)|) (28) 

In the presence of a source of wave energy, continu- 
ous shearing/straining experienced by wave packets along 
their trajectory leads to a cascade of their energy, and 
hence, of the energy of the global wave perturbation. This 
cascade can be studied in two different ways. In the first, 
there is no steady source, but only an initial distribution 
of wave packets and energy in space. For this initial value 
problem, the interest is in the subsequent transient evo- 
lution of the wave energy spectrum Fi{k 1 t) with initial 
condition F/(fc, t = 0) =J2 e j( t = 0)5(k- \k 3 (t = 0)|). As 
time increases, each wave number kj increases in time (on 
average), evolves toward the dissipation range and, there- 
fore, the wave energy decays. A second approach consists 
in taking into account continuous injection of wave energy 
given by the ongoing perturbations. Therefore, wave pack- 
ets and energy are supposed to be continually injected at 
intermediate k values, around ko, a typical injection scale, 
with L^ 1 <C fco <C kd and with kd the very small dis- 
sipative cutoff and L the global length scale associated 
with the magnetic field. In this case, electromechanical 
energy injection, say the photospheric drive, becomes bal- 
anced by heating and dissipation at high k, and a time- 
averaged stationary spectrum F$(k) builds- up. We turn to 
the study of the wave energy spectrum that results from 
quasi-uniform Alfvenic strain experienced by wave-packets 
as they propagate in a field with destroyed magnetic sur- 
faces, i.e. along chaotic magnetic field lines. This config- 
uration, which is a paradigm in transport theory, for the 
calculation of anomalous diffusion coefficients (Rechester 
and Rosenbluth 1977), was first considered in the con- 
text of coronal wave dissipation by Similon and Sudan 
(1989), although they were only interested in the mixing 
time scale. Calling e the input wave energy flux, dimen- 
sional analysis yields 

e oc V A kF s (k), (29) 

with Fs(k) the stationary wave energy spectrum. Hence, 

F s (k) cx k~ l V A l e, (30) 

in the range where dissipation is negligible, i.e. for fc < 
kd = {V A /n)^ 2 . This relation §U$i for the wave-energy, 
which is again derived below from a kinetic equation, is 
the same as the Batchelor's spectrum (Batchelor 1959) 
for scalar variance, obtained for a flow with uniformly dis- 
tributed rate of strain, i.e. with chaotic streamlines. 
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An other derivation of ([5U|) follows from writing down 
the evolution equation for the wave energy spectrum 
F(k,t). Taking the time derivative of Eq. (f2"5)) . we obtain 



= -r,k 2 F(k,t) + ^-[V A kF(k,t)], 



(31) 



The first term on the right-hand side represents the effect 
of dissipation, e = —rjk 2 e, and vanishes for n = 0. The 
second term describes transport in k-space for wave pack- 
ets that experience uniform rate of strain, V' A — const, 
i.e. k(t) — V A k. The stationary solution of (|3Tj) for 77 = 
is again F s (k) oc k^ 1 . For 77 ^ 0, the stationary solution 
to Eg. pip , with energy source e added at wave numbers 
where wave energy is fed into the system, is 



F s (k) 



V A k 



exp(~nk 2 1 a) 



(32) 



This shows the existence of a gaussian dissipative fall-off, 
for F s (k), at large k > kd- 

Let us note that while the derivation of the kinetic 
equation (|3Tj) is based on the WKB approximation, the 
dimensional analysis obtention of the k~ l stationary spec- 
trum does not rely on the assumption of a separation 
of scales. The latter is based only on the fact that the 
wave energy transfer time is a constant, independent of 
the scale, given by V A , the quasi-uniform rate of strain 
experience by the waves. 

Finally, we go back to the initial value problem, for 
transient forcing. Discarding dissipation in Eq. ([3"Tj) . and 
changing to the new variables, £ = ln(fc) and G — kF(k), 
we can write Eq. (|3"T)) as 



dG 
~dt 



= V: 



,dG 



(33) 



Therefore, this transport equation reveals that, if an in- 
put wave energy pulse, Fj(k,t = 0), is initially centered 
at intermediate k values, say k — kg, it is then carried 
conservatively at constant speed a = V A , on a logarith- 
mic scale, until it reaches the dissipative range k — kd- 
Only there it experiences fast damping, and finally con- 
verts its energy into heat. This spectral viewpoint also 
contains the time scale involved before dissipation occurs 
: it is the travel time in the energy conserving range in k- 
space, r = {\nk d -\nk Q )/V A = (2/V A )ln{V A /nk 2 ) oc In S, 

Which is T st rain- 

Before concluding, few comments are due. First, let 
us note that the magnitude of the energy flux e depends 
on the amplitude of the drive, but not on its frequency, 
see e.g. Velli (2003). However, the location of the energy 
source in fc-space depends on the range of driving frequen- 
cies. Recall that since the dynamical equations are linear, 
consideration of a broad frequency drive does not alter 
the shape of the spectrum in the gap between the largest 
input wavenumber and the dissipative one. The very ex- 
istence of such conservative range with neither source nor 
dissipation is suggested by the following reasoning. Since 
the wave-number associated with the dissipative scale is 
very large, kd = (YaM) 1 ^ 1; the condition that the 



latter is also much larger than the largest input wave- 
number kn max = uj/Va, i-e. k^ max <c kd, is equivalent to 
lot a < L(V A /n) 1/2 ~ S 1/2 . This condition can easily hold 
for also for waves with ujt a 3> 1 or equivalently, fc|| 3> L^ 1 , 
for which the WKB approximation holds. In their origi- 
nal work, Similon and Sudan (1989), were able to relax 
the condition ku 3> making it possible to study the 
behavior of the lower frequency fluctuations and "non- 
WKB" effects. Their tactic is entirely analogous to the 
reduced magnetohydrodynamics (RMHD) approximation, 
although they refer to the so-called "ballooning represen- 
tation" in fusion studies. By doing so, they only restrict 
the eikonal representation to the perpendicular direction. 
Low frequency motion at the boundary of a unidirectional 
magnetic field is well known to produce turbulence and 
energy cascade leading to a power-law energy spectrum 
[see for instance two recent works, Rapazzo et al (2007), 
Dmitruk et al. (2003) and references therein]. For realis- 
tic forcing, i.e. small amplitudes of the Alfven waves that 
are launched into the corona, 5u ~ 10~ 3 Va (Acton et al. 
1981), numerically obtained spectra are generally found 
to be quite steep. The reason resides in the weakness of 
the nonlinear interactions. However, since the background 
large scale field is taken as a constant, only the nonlin- 
earity can there provide the necessary transfer of energy 
to small dissipative scales. The linear mixing process with 
its associated energy flux is likely to alter energy trans- 
fer rates and hence spectra in inhomogeneous magnetic 
backgrounds. This is suggested by the numerical studies 
by Malara et al. (1992), Einaudi et al. (1996), who showed 
that scattering of Alfven waves with the background inho- 
mogeneity, for a population of waves initially propagating 
say, upward, can act as a source for downward-propagating 
waves, and hence as a trigger for the occurrence of non- 
linear interactions. We did not consider such effect in the 
present work and we are not aware of any MHD turbu- 
lence theory that accounts for how the linear and non- 
linear mechanisms of energy flux combine together. We 
find this subject an interesting area for further research. 



We finally mention that the mixing effect considered 
in the present article, is the root of a proposed geomet- 
rical interpretation for the consequence of non-linear in- 
teractions among counter-propagating wave packets. It 
is clearly illustrated by Fig.l in Maron and Goldrcich 
(2001). Indeed, take a sample of originally straight ver- 
tical field lines which are perturbed by, say, downward- 
propagating waves. As a consequence, the field lines wan- 
der around each other. Therefore, if a localized wave 
packet is launched and propagates upward in this braided 
field, it will be distorted [also compare with Fig. 2 in 
Similon and Sudan (1989) or the figure in Rechester and 
Rosenbluth (1977)]. After several of these collisions, dis- 
sipation takes over and the wave packet mixes in its envi- 
ronment. 
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5. Conclusions 

In this work, we consider the wave energy flux to small 
scales which results from the distortion of Alfven waves 
propagating in an inhomogeneous mean magnetic field Bq. 
This problem is related to turbulence because the exis- 
tence of a mean field may provide separation between lin- 
ear and non-linear time-scales. In an inhomogeneous mag- 
netic field, Alfven wave packets disperse and are stretched 
by an Alfvenic flow = Bo/(^o/ 3 ) 1/ ' 2 , whose streamlines 
are also the magnetics field lines. In a plasma with uni- 
form density p, this flow is incompressible, i.e. V.V^ = 0, 
and its stretching properties are thus only related to mag- 
netic field lines dispersion. While the wave packets trans- 
port their energy at the group velocity V^, and hence 
while they disperse, in physical space, as the field lines 
do, they also experience continuous shearing and strain- 
ing, i.e. they tend to mix. This gives rise to a cascade of 
the wave energy in /c-space. As is well known for mixing, 
the efficiency of stretching relies on the local exponential 
separation of nearby trajectories (Ottino 1988). Since for 
a continuous external wave energy drive, the cascade of 
this energy ultimately leads to a balance between source 
and heating in the dissipative range, a stationary wave 
energy spectrum can build up even in the absence of sig- 
nificant nonlinear interactions among the waves. A phe- 
nomenon which is absent in a constant mean magnetic 
field. Assuming the existence of a conservative range, we 
determined the wave-energy spectrum for the special case 
of wave packets propagating in a " fully" chaotic magnetic 
field. It means that the distribution of "Lyapunov expo- 
nents" , representing the rate of local exponentiation of 
nearby field lines, and hence the degree of stretching of 
wave packets, is strongly peaked at an average value a. 
This situation, which was first considered by Similon and 
Sudan (1989), for dissipation of Alfven waves corresponds 
to chaotic mixing by flows with uniformly distributed rate 
of strain (Batchelor 1959, Antonsen et al. 1996). 
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